Recovering MOND from extended metric theories of gravity 



O 

(N 
-t— > 
O 

O 



T. Bernal 1 ^ S. Capozziello 2 ' 3 Q J.C. Hidalgo 1 ' 4 ! and S. Mendoza 1 ! 

1 InsUtuto de Astronomia, Universidad National Autdnoma de Mexico, AP 70-264, Distrito Federal 04510, Mexico 
2 Dipartirnento di Scienze Fisiche, Universita degli Studi di Napoli "Federico II", 
Complesso Universitario di Monte Sant'Angelo, Edificio N, via Cinthia, 80126, Napoli, Italy 
3 1. N.F.N. - Sezione di Napoli, Complesso Universitario di Monte Sant'Angelo, Edificio G, via Cinthia, 80126, Napoli, Itah 
4 Departamento de Fisica, Instituto Nacional de Investigaciones Nucleares, AP 18-1027, Distrito Federal 11801, Mexico 

(Dated: October 28, 2011) 

We show that the Modified Newtonian Dynamics (MOND) regime can be fully recovered as the 
weak-field limit of a particular theory of gravity formulated in the metric approach. This is possible 
when Milgrom's acceleration constant is taken as a fundamental quantity which couples to the theory 
in a very consistent manner. As a consequence, the scale invariance of the gravitational interaction 
is naturally broken. In this sense, Newtonian gravity is the weak-field limit of general relativity 
and MOND is the weak-field limit of that particular extended theory of gravity. We also prove 
that a Noether's symmetry approach to the problem yields a conserved quantity coherent with this 
relativistic MONDian extension. 
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I. INTRODUCTION 

Milgrom [H developed a non-relativistic theory of 
gravity in order to explain observed flat rotation curves 
of spiral galaxies. This Modified Newtonian Dynamics 
(MOND) theory has proved useful in explaining a great 
variety of astronomical phenomena without requiring the 
presence of a dark matter component (see e.g. [H, H| and 
references therein). 

As explained by an extended Newtonian approach to 
gravity, Mendoza et al. [5( showed that the key feature of 
MOND is the introduction of a fundamental acceleration 
ao in the theory, which on itself at the non-relativistic 
level, makes gravity scale dependent. 

Through the years, finding the relativistic extension 
of MOND has become a big challenge. The most suc- 
cessful attempt was proposed by Bekenstein Q who for- 
mulated a Tensor- Vector-Scalar (TeVeS) relativistic the- 
ory of MOND. This approach requires tensor, vector and 
scalar fields to achieve a self-consistent description. How- 
ever, the many cumbersome mathematical complications 
that TeVeS present are evident. Furthermore, it cannot 
reproduce crucial astrophysical phenomena (see e.g. 0). 

On the other hand, since extended metric f(R) the- 
ories of gravity have proved very successful on a wide 
variety of cosmological scenarios [see e.g. [8l-[TTl and ref- 
erences therein] it is natural to seek a relativistic gener- 
alisation of MOND in this direction, which has not de- 
served too much attention since Soussa [l2[ and Soussa 
and Woddard [l3[ developed a no-go theorem which pre- 
vented all metric f(R) theories of gravity to become rel- 
ativistic candidates for MOND. However, Mendoza and 
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Rosas-Guevara [14[ showed counterexamples of this no- 
go theorem disproving its general validity. Furthermore, 
the works by Capozziello et al. fl5l - fll| and Sobouti [l8| 
made clear that particular f(R) models are capable of 
explaining phenomena usually ascribed to MOND. Up 
to now, the challenge has been to formulate a suitable 
f(R) theory able to converge to standard MOND in the 
non-relativistic regime. 

In this article, we show how a particular f(R) metric 
theory of gravity, derived from first principles, is capable 
of reproducing MOND when its non-relativistic regime is 
reached. To do so, in Section [XT] we set the foundations 
of a metric theory of gravity with the use of correct di- 
mensional quantities. Then, in Section IIIII we solve the 
problem of a point mass source on a static space-time and 
find, at first order of approximation, the particular form 
of the function f(R). In Section llVl the integration con- 
stants of the theory are fixed by solving the same problem 
in the formalism of metric perturbations. In SectionlVlwe 
show that the existence of a Noether symmetry confirms 
the appearance of a characteristic scale of the problem in 
the MONDian regime. Finally in Section fVTl wc comment 
on the obtained results. 



II. DIMENSIONAL GROUNDING 

Let us assume that a point mass M located at the 
origin of coordinates generates a relativistic gravitational 
field in the MONDian regime and that a metric formalism 
describes the field equations. 

This problem is characterised by the following quanti- 
ties: the speed of light in vacuum c, the mass M of the 
central object generating the gravitational field, New- 
ton's constant of gravity G and Milgrom's acceleration 
constant oq. With these parameters, two "fundamental 
lengths" can be built: 
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The gravitational radius r g is a length that appears once 
rclativistic effects are introduced on a theory of grav- 
ity. The mass- length scale l M , as described in [5], is 
a characteristic length which appears on a gravitational 
theory when MONDian effects are to be taken into ac- 
count (for consistency we note here that a third length, 
A := l]rf/r g = c 2 /ao, that does not contain the mass, can 
be constructed as a combination of the previous two). 

As discussed by Mendoza et al. 5], in the non- 
rclativistic regime, a test particle located at the radial 
coordinate r from the origin will obey the MONDian dy- 
namics when Im/t 1- When Iw/r 3> 1 the gravi- 
tational field is Newtonian. As such, when relativistic 
effects are taken into account for the gravitational field, 
then standard general relativity should be recovered in 
the limit Im/t 1, and a relativistic version of MOND 
should be obtained when Im/t "C 1. This shows that 
the pursue of a complete metric description leads one to 
consider the scale-dependence of gravity. 

The length scales presented in equation (JTJ must some- 
how appear in a relativistic theory of gravity which ac- 
cepts the fundamental nature of the constant a®. For 
example, in the metric formalism, a generalised Hilbert 
action Sh can be written in the following way: 
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which slightly differs from its traditional form (see e.g. 

iii i3) 



(3) 



since we have introduced the following dimensionless 
quantity: 



X 



L M R, 



(4) 



where R is Ricci's scalar and Lm defines a length fixed 
by the parameters of the theory. The explicit form of the 
length Lm has to be obtained once a certain known limit 
of the theory is taken, usually a non-relativistic limit. 
Note that the definition of x gives a correct dimensional 
character to the action (J2J), something that is not com- 
pletely clear in all previous works dealing with a metric 
description of the gravitational field. For f{\) — X the 
standard Einstein-Hubert action is obtained. 

On the other hand, the matter action has its usual 
form, 



Sm ~ 2~c 



-gd x, 



(5) 



with £ m the Lagrangian density of the system. The null 
variations of the complete action, i.e. S (Sh + S m ) = 0, 
yield the following field equations: 
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where the dimensionless Ricci tensor \y.v is given by: 



(7) 



and R^ is the standard Ricci tensor. The Laplace- 
Beltrami operator has been written as A := V Q V a 
and the prime denotes derivative with respect to its 
argument. The energy- momentum tensor T^„ is de- 
fined through the following standard relation: 5S m = 
— (l/2c) T a p 8g a P . In here and in what follows, we choose 
a (+,—,—,—) signature for the metric g^ v and use Ein- 
stein's summation convention over repeated indices. 
The trace of equation ^ is: 



/'(x)x-2/( X ) + 3L 2 M A/'(x) = 



8irGLl 



M 



T, (8) 



where T := T°. 

Since we are only interested on the gravitational field 
produced by a point mass source located at the origin, 
then the mass density p is given by 



MS(r), 



(9) 



where 5(r) represents the three dimensional Dirac delta 
function. With this, it follows that the only non-zero 
component of the energy-momentum tensor is given by 



pc 



l M8{r) 



(10) 



A point mass distribution generates a stationary spher- 
ically symmetric space-time and so, the trace equa- 
tion (JSJ contains all the relevant information relating the 
field equations. In what follows we assume a power law 
form for the function f(x), i-c. 



fix) = x b 



(11) 



III. ORDER OF MAGNITUDE APPROACH 

Let us first analyse the problem described in the previ- 
ous Section by performing an order of magnitude approx- 
imation of the trace equation (|8[). Under these circum- 
stances, d/dx ~ 1/Xj A ~ — 1/r 2 and the mass density 
p m M/r 3 . This approximation implies that equation ([5]) 
takes the following form: 
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Note that the second term on the left-hand side of 
equation (1121) is much greater than the first term when 
the following condition is satisfied: 



Rr 2 < 
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At the same order of approximation, Ricci's scalar R w 
k = R~ 2 , where k is the Gaussian curvature of space and 
R c its radius of curvature and so, relation (|13[) essentially 
means that 



R c > r. 



(14) 



In other words, the second term on the left-hand side 
of equation (fT2"j) dominates the first one when the lo- 
cal radius of curvature of space is much grater than the 
characteristic length r. This should occur in the weak- 
field regime, where MONDian effects are expected. For a 
metric description of gravity, this limit must correspond 
to the relativistic regime of MOND. In this article we 
will only deal with this approximation. At the end of 
the current Section we show an equivalent relation to in- 
equality (|14[) which has a more physical meaning. 

Under assumption ([TO , equation (fl~2"|) takes the follow- 
ing form: 



This last equation converges to a MOND-like acceler- 
ation acxl/rif& — 2 = — (6— 1), i.e. when 



b = 3/2. 



(19) 



Also, at the lowest order of approximation, in the extreme 
non-relativistic limit, the velocity of light c should not 
appear on equation (|18|) and so, the only possibility is 
that Lm depends on a power of c, i.e. 



L 



M 



(X c 



(4-2b)/(6-l) 



c 2 , 



and so, 



Lm oc c . 



(20) 

As discussed in Section UH the length Lm must be con- 
structed by fundamental parameters describing the the- 
ory of gravity and so, let us assume that 



L M = W M , 



with 



a + P = 1, 



(21) 



where the constant of proportionality £ is a dimensionless 
number of order one that will be formally obtained in 
Section IIV1 Substituting equation (f2"Tj) and the value 
obtained in (IT91 into relation ([2"U)l it then follows that 



= 0=1/2, 



i.e. 



J->M ~ r i M 



(22) 



If we now substitute this last result and relation (|19p 
in equation (IT8l) it follows that 



8ttGA/ 



(15) 



We now recall the well known relation followed by the 
Ricci scalar at second order of approximation at the non- 
relativistic level 12011: 



R = -- 

c z 



-V-a, 

c z 



(16) 



where the negative gradients of the gravitational poten- 
tial <j) provide the acceleration a := — V0 felt by a test 
particle on a non-relativistic gravitational field. At order 
of magnitude, equation (fTU)) can be approximated as 



(a a GM) 1/2 



(23) 



which is the traditional form of MOND in spherical sym- 
metry (see e.g. [1, 0, Hl| and references therein). Also, 
the results of equation (|23|) in (fTT|) mean that 



R 



Im r 2 ' 



(24) 



and so, inequality (|14[) is equivalent to 



R 



2(b 



2a 



(17) 



Substitution of this last equation on relation (I15|) gives 



2 r {8ttGM\ 



1/(6-1) 



err 

~2L~% { 3bc 2 r ) 

_ c (26-4)/(6-l) r (6-2)/(6-l) i -2 ^^1/(6-1) 



(18) 



Im > r g 



(25) 



The regime imposed by equation (|25p is precisely the one 
for which MONDian effects should appear in a relativistic 
theory of gravity. This is an expected generalisation of 
the results presented by Mendoza et al. [|| in the weak 
field limit regime for which Im *C r and so, combining 
this with equation ([231) yields r ^> Im ^> r g . In this 
connection, we also note that Newton's theory of gravity 
is recovered in the limit I m 3> r 3> r„ . 
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IV. WEAK FIELD LIMIT APPROACH 

We now use the trace ([8]) to the lowest order of per- 
turbation. Results of this perturbation in the Newto- 
nian limit for other metric theories of gravity have been 
reported by [22|, However, since we are interested 
in the lowest order of approximation in the MONDian 
regime, we expect different results. 

Under the assumption of spherical symmetry for a 
static space-time, with a diagonal metric given by 



The trivial solution of the previous Laplace-like equation 
occurs when the argument of the Laplacian is equal to 
zero, and so 

f'Od = -^ N - ( 31 ) 

Substitution of equations ((11]), CEH]), (|2"T]). and ^ in 
relation (f3~Tj) leads to 



2(j) 

900 = 1 + -o- , 

922 = ~r 2 , 



then 



9u 

933 



-r 2 sin 2 I 



(26) 



R 



M 



(32) 



where we have used the fact that for a point mass M, 
the Newtonian potential is given by: 



A/'(%) = 



1 



d 



—g dx^ 
1 d 



di 



v-gg 
2 df(x) 

Or 



dx v 



(27) 



-v 2 /'(x), 



at the lowest order of perturbation in (j>/c 2 . In this case, 
when condition (TI4")) or equivalently equation (|2~51 is sat- 
isfied, the trace © of the field equations is given by: 



3V 2 /'( X ) 



8irG 



(28) 



Note that this equation can also be obtained by perform- 
ing a direct perturbation of ®. This is so because in 
the weak field limit, the field equations (JXJ) are studied 
at orders of powers of c -2 [13, HH ■ At the lowest zeroth- 
order of the perturbation, the Ricci scalar = every- 
where and so, it describes a flat space-time. At the next 
second perturbation order O(x) = O(-^m) + C( (2) i?), 
where ^R represents Ricci's scalar at the second order 
of the perturbation. Since 0(Lm) = 1 according to re- 
lation (|2"0")) . then 0(x) = 4. If we now use the power 
law relation (fTTj) with (|19p and perturb the trace (|S]), it 
follows that the first two terms on the right-hand side 
are of order 0(x 3 ^ 2 ) = 6, while the remaining two terms 
are of order 0(L 2 M x x ^ 2 ) = 4 and so, it follows that the 
trace (|5]) at its lowest non-zero perturbation order is ex- 
actly relation ([28)) . 

Equation (f2"5|) can be integrated straightforward using 
the standard Poisson equation of Newtonian gravity: 



G- 



(33) 



We now substitute relation (ITBl into equation (l3"2"j) to 
obtain: 



(34) 



Integrating this equation over a spherical volume of 
radius r and using Gauss's theorem on the left-hand side, 
we obtain 



'2V2 N 
9C , 



c 2 r„\ 1 



hi J r ' 



(35) 



and so, by choosing 



2y/2 



we reach the MOND acceleration limit 



(36) 



(a GM) 1/2 



(37) 



V 2 </> N = AnGp, (29) 

for the Newtonian potential fa- Substitution of this re- 
lation on (|2"8"|) yields 

V 2 (f'(x) + ^n) = 0. (30) 



at the lowest perturbation order of the theory. 

The formulation of the theory described so far implies 
that, as soon as we relax the hypothesis for which the 
gravitational action is strictly the one described by stan- 
dard general relativity, new characteristic lengths have 
to be considered. In the following Section, we show that 
such scales coherently appear as constants of motion of 
the problem. 
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NOETHER'S SYMMETRIES 



The above results are based on the fact that we are as- 
suming the power law relation (ITT1) for the gravitational 
action. In particular, for f(R) oc i? 3 / 2 the MOND ac- 
celeration regime is recovered. This particular theory 
also admits exact cosmological solutions where a matter 
dominated era evolves towards the accelerated universe 
observed today [24| and recovers the observed dynamics 
at astrophysical scales (25|. Also, this solution is partic- 
ularly relevant since its conformal transformation is ex- 
actly invertible, as shown in [26]. Noether's symmetries 
give rise to conserved quantities that are directly related 
to characteristic length scales 



T 2 

M 



24] 



In order to develop a Noether's approach to the prob- 
lem, let us consider /(%) gravity in static spherical sym- 
metry, following the same ideas as the ones exposed in 
Section lTVl The spherical point-like /(x) Lagrangian can 
be obtained by imposing the spherical symmetry directly 
into the action J5}. As a consequence, the infinite num- 
ber of degrees of freedom of the original held theory will 
be reduced to a finite number. The technique is based 
on the choice of a suitable Lagrange multiplier defined 
by assuming the known explicit form of Ricci's scalar R 

m 

The static spherically symmetric metric can be ex- 
pressed as 



ds 2 = A(r)c 2 dt 2 - B(r)dr 2 - C(r)dft 2 , (38) 

where dil 2 := d# 2 +sm 2 9dip 2 is the angular displacement 
and C(r) :— r 2 . Since we arc interested in the weak- 
field limit approach of the /(x) theory, let us assume 
B(r) = 1. 

The point-like f(x) Lagrangian L is obtained by 
rewriting the action @ as 



S = 



l6irGL 2 M 



[f(x)~Hx-x)}V^gd 4 x, (39) 



where A is a Lagrange multiplier and \ = L 2 M R, for 
the known Ricci scalar R expressed in terms of the met- 
ric (HHJ) with B(r) = 1: 



R- A " 



2C" 



A'C 



A' 2 C' 2 2 
2A 2 2C 2 C [ j 



C AC 

In the previous equation, the prime denotes the derivative 
with respect to the radial coordinate r. Variations of the 
action with respect to \ gi ye the explicit form of the 
Lagrange multiplier: 

d X Jx 

Substituting this result in the action (|39[) and elimi- 
nating the boundary terms (cf. [27]), the point -like La- 
grangian is obtained: 



2Af xx C' X ' 



-VA[(2L 2 M + C X )f x -Cf] 



(41) 



Note that this Lagrangian is canonical since only the gen- 
eralised positions, q — (A,C,x), and their generalised 
velocities q' — (A',C',x')> appear explicitly. 

In the MONDian regime, where equation (| X4[) holds, 
the last two terms on the right-hand side of equation (|4Tj) 
arc of order Cf ~ Cx 3 ^ 2 and so, both are much smaller 
than L 2 M f x . This statement is also true at the lowest 
non-zero order of perturbation, since all terms on the 
right-hand side of equation (|4"Tj) are of order 4, except 
the last two which are of order 6. With this, it follows 
that in the MONDian regime the Lagrangian (l4"Tj) can be 
written as 



T 2 



Ah 



Va 

2Af xx C' X ' 



2C 



C' 2 + f x A'C + Cf xx A' x '- 



2A 



(42) 



We now search for symmetries related to the cyclic 
variables and then reduce the dynamics. According to 
Noether's theorem, the existence of symmetry properties 
for the dynamics described by the Lagr ang ian implies 
the existence of conserved quantities [28l43C| . In princi- 
ple, this approach allows to select particular f(x) gravity 
models compatible with spherical symmetry. 

A conserved quantity exists if the Lie derivative of the 
Lagrangian (14"2")) along the vector field X vanishes: 



£xL = 



*'V„<L = 0, 



(43) 



for 2 = 1,2, 3, in the configuration space (A, C, x)- Solv- 
ing equation (1431) means to find out the functions on 
which constitute the Noether vector [2^, (29[. However, 
the relation (|43|) implicitly depends on the form of f(x) 
and then, by solving it, we also get f(x) models com- 
patible with spherical symmetry. On the other hand, by 
fixing the f(x) function, we can explicitly solve (|43|) . In 
principle, the same procedure can be worked out any time 
Noether's symmetries are identified [3lj . 

The general form of the Noether vector is given by the 
solution of the Killing equations for the components oti 
of the vector a in a flat space-time (cf. [321]): 



ax — k x A +pi, 
a 2 = k 2 C + P2, 
"3 = hx+P3, 



(44) 



with hi and pi constants. Using the power law (jlll) and 
the general solution (|4"4")l in equation (|4"3"]l . it follows that 
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a = \ 2(1 - b)kA, 0, k X j , (45) 

where k is an integration constant and we have assumed 
in the calculations that 6 ^ 1, yf ^ and 2AC'+A'C ^ 
in order to obtain solutions different from general relativ- 
ity (cf. Section |TTT|. For this case, the related constant 
of motion, Eo, is given by 

So = QjjV^L, 

= L 2 M 6(6 - l)kA-V 2 C X h - 2 [2(6 - 1)A X < - A'x] . 

(46) 

In general relativity, where 6=1, Noether's sym- 
metry approach gives Eo = 2r g , which is exactly the 
Schwarzschild radius (cf. [13]). On the other hand, in 
the MONDian regime, where equation (14T>1) is valid for 
6 = 3/2, C(r) = r 2 and at the lowest order of perturba- 
tion A(r) = 1 + 2(p/c 2 , the related constant of motion is 
given by 

E = \krll M . (47) 

The key point in this relation is that it includes two 
characteristic lengths, namely the mass length-scale Im 
and the gravitational radius r g . This is coherent with 
the results discussed on the previous Sections, since both 
characteristic lengths must appear on a correct relativis- 
tic metric theory of MOND. 

VI. DISCUSSION 

We have shown that a metric theory of gravity f(x) — 
X 3 / 2 , where the dimensionless Ricci scalar X is given by 
equation (|4|, converges to MOND in the non-relativistic 
regime. We note that a previous attempt in this di- 
rection was made in the B.Sc. thesis of Rosas-Guevara 



[33} . where a metric f(R) oc i? 3 / 2 theory was proposed 
to account for relativistic effects of MOND. However, 
the correct approximation (f25j) was never introduced on 
that analysis and therefore, the MONDian limit was not 
achieved on that work. 

We have also shown that the appearance of a new char- 
acteristic length related to MOND's acceleration is coher- 
ent with Noether's symmetries related to the problem. 

The metric theory of gravity presented here is by no 
means a complete description at all scales of gravitation. 
It only deals with the MONDian regime of gravity, i.e. 
when equation (|25[) is valid. In other words, our de- 
scription breaks the scale invariance of gravity in a more 
general way than the one described in |5| . 

The mass dependence of x means that the mass needs 
to appear on Hilbert's action ©. This is traditionally 
not the case, since that action is thought to be purely a 
function of the geometry of space-time due to the pres- 
ence of mass and energy. However, it was Sobouti [l8[ 
who first encountered this peculiarity in the Hilbert ac- 
tion when dealing with a metric generalisation of MOND. 
Following his remarks [l8j . one should not be surprised 
if some of the commonly accepted notions, even at the 
fundamental level of the action, require generalisations 
and re-thinking. An extended metric theory of gravity 
goes beyond the traditional general relativity ideas and 
in this way, we probably need to change our standard 
view of its fundamental principles. 
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